Theory of small charge solitons in one-dimensional arrays of Josephson junctions 



ON 

o 
o 

(N 
> 

o 



c ■ 
o ■ 

o , 

I 

=3 : 

c3 ■ 



O 

i o i 

(N 
> 

in 

On 

(N 

l> 
O 
On 
O 



X 



Stephan Rachel and Alexander Shnirman 
Institut fur Theorie der Kondensierten Materie and DFG Center for Functional Nanostructures, 
Karlsruhe Institute of Technology, D-76128 Karlsruhe, Germany 

We identify and investigate the new parameter regime of small charge solitons in one-dimensional 
arrays of Josephson junctions. We obtain the dispersion relation of the soliton and show that it 
unexpectedly flattens in the outer region of the Brillouin zone. We demonstrate Lorentz contraction 
of the soliton in the middle of the Brillouin zone as well as broadening of the soliton in the flat band 
regime. 
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Charge solitons in one-dimensional (ID) arrays of tun- 
nel junctions in the Coulomb blockade regime were in- 
troduced about twenty years ago [l|, 0] and are being 
studied ever since (see, e.g., Ref. Q). Hermon et al. 0| 
studied a onc-dimcnsional array of Josephson junctions 
(JJs). It was shown that, if the grains have a large ki- 
netic (or geometric) inductance, the system's dynamics 
are governed by the sine-Gordon model and, therefore, 
kink- like topological excitations, i.e., charge solitons, are 
the charge carriers. Simultaneous experiments by Hav- 
iland and Delsing [f| demonstrated the Coulomb block- 
ade in ID arrays of JJs consistent with the existence of 
charge solitons. In the later experiments of Haviland's 
group 0, 0] considerable hysteresis in the I-V character- 
istic of the array was observed and attributed to a very 
large kinetic inductance. The physical origin of this in- 
ductance remained unclear. A few years later, Zorin [f| 
pointed out that a current biased small-capacitance JJ 
develops an inductive response on top of the capacitive 
one. This phenomenon was called Block inductance. A 
closely related inductive coupling between two charge 
qubits was studied in Ref. . It is still not clear if Bloch 
inductance could support the dynamics of charge soli- 
tons. 

In this paper, we identify the new parameter regime 
within the Coulomb blockade (insulating) phase of a ID 
array of coupled JJs. It is defined by the condition 
AEj > Ec > Ej, where Ec and Ej are the charging 
and the Josephson energies of the junction, respcctivly, 
and A is the bare screening length (measured in number 
of junctions). In this regime we investigate the dynamics 
of charge solitons and demonstrate two surprising fea- 
tures: i) flattening of the dispersion relation in the outer 
region of the Brillouin zone; ii) broadening of the soli- 
ton in the flat band regime in contrast to the expected 
and observed Lorenz contraction in the regime of regular 
dispersion relation. We believe these results might open 
the wav to the explanation of the experimental data of 
Refs. " 
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The paper is organized as follows. In order to shed 
light on the previous studies of charge solitons in terms 
of the relativistic sine-Gordon equation and to facilitate 
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FIG. 1: Josephson junction array. 

the interpretation of our new results we, first, formulate 
the mean-field approach. Then we develop a many-body 
tight binding technique which leads to the new results. 

The system considered is shown in Fig.[TJ The grains 
are connected by JJs of capacitance C (typically 1 fF) 
and each grain has a capacitance Co to the ground (typ- 
ically 5 — 20 aF). The kinetic or geometric inductance of 
the grains Lq is included to simplify the mean-field treat- 
ment but it is later assumed to be vanishingly small. We 
derive the following Hamiltonian: 
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Here m r is the number of Cooper pairs that have 
tunneled through junction number r. The continuous 
polarization charge Q r = J2r'<r 1^° corresponds to 
the integral of current flown into junction number r. 
The commutation relations read [$ r ,Qr'] = ift5 rr > and 
[m r ,e^'] = e^Sry. 

Mean field approach. In the mean field approxima- 
tion we treat the dynamical variables Q r as c-numbers, 
Qr — * (Qr(t))- The Hartree-like wave function can be 
written as a product of single junction states, \I/({m}) = 
Tlr *Q r ( m r)- Here 'I'q(to) is the (ground) state of a sin- 
gle junction with Hamiltonian 
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The self-consistency condition is derived by averaging the 
equation of motion for the variables Q r : 

2Qr Qj.J r -\ Qr — 1 
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where V r = (Q r — 2em r )/C is the average voltage on 
the junction r. For static Q r and at zero temperature 
V r = dE (Q r )/dQ r , where E (Q) is the lowest energy 
band of Hamiltonian . Zorin [8[ derived an additional 
inductive contribution to the voltage on the junction: 
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LB(Qr)Qr, where Lb(Q) is the Bloch 



inductance. Then, Eq. ([3]) reads 
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where L c g = Lq + LsiQr)- We observe that the induc- 
tance Lq is superseded by the Bloch inductance and we 
can safely assume Lq = 0. 

For the case of Q-hidependent inductance L e ff , Eq. (|4|) 
was studied in Ref. Q (there it was assumed that the 
inductance is dominated by the kinetic inductance of 
the superconducting islands). Eq. (g]) is, then, a dis- 
crete analog of the relativistic sine-Gordon equation and 
it possesses topological solitons which describe the prop- 
agation of Cooper pairs through the array. As usual in 
relativistic physics, a soliton is subject to the Lorentz 
contraction, i.e., its length reduces as its velocity grows 
(see Ref. [H and references therein). 

Investigation of the case of Q-dcpendent inductance 
(the Bloch inductance is a rapidly varying function of Q 
in the regime Ec > Ej) is still pending. We just note 
here that one could expect 0] the effective Lagrangian 
of a Q-biascd Joscphson junction to have the form C = 
(1/2)Lb(Q)Q 2 —Eq(Q). Thenthevoltageonthejunction 
r would be given by V r = fg& + L B (Q r )Q r + \^Ql- 
Thus, Eq. (j4|) might need to be further modified. In this 
paper we do not pursue further the mean field analy- 
sis but rather concentrate on an alternative approach of 
tight binding treatment of various charge configurations. 

Charge configurations. For Lq — > the polariza- 
tion charges Q r are enslaved to the discrete charges 
m r (the charge that have tunneled through junction 
r). If the charge configuration {m r } is given, then the 
polarization charges {Q r } are found from Q^~ 2em ^ 



2Q r 



C 



— ''c„~ 1 = 0. Equivalently one can consider is- 
land charges n r = m r — m r+1 and obtain the charging 
energy of the array (see, e.g., Ref. (loj ) 



H c = - ^ U(r - r') n r n r 
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where A = \J C/Cq is the screening length and Ec = 
(2e) 2 /(2C) is the charging energy of a single junction. 
The Joscphson term in the Hamiltonian connects the 
charge configurations which differ by one Cooper pair 



being transported through one junction. For A ^> 1, the 
charging energy reads U{r) w AEc exp (— |r|/A). 

Charge states nomenclature. We consider the sector 
of the Hilbert space with exactly one extra Cooper pair 
in the array, i.e., ^2 r n r = 1. The simplest representa- 
tive of this sector is the state in which the extra Cooper 
pair resides on island R and all other islands are neutral. 
We denote this state |. . . 1 R . . .) = \R). The charg- 
ing energy of \R) is given by \U{0) = E w AE c /2. 
This is a rather high energy, in case of Cq — > it 
is in fact infinite (proportional to the system size [lit), 
and this is approximately the energy one has to in- 
vest in order to insert the Cooper pair into the ar- 
ray. There exists, however, other charge configurations 
in the single Cooper pair sector, i.e., the ones with 
charge-anti-charge pairs induced in the vicinity of the 
first Cooper pair. The first example is the configura- 
tion [...001 - l fl 100...) = |-R;1,1), where charge -1 
resides on island R while charges +1 reside on the neigh- 
boring islands R—l and R+l. Its charging energy is given 
by E a + E 1A , where E x . x = U(0)-2U(1) + U(2) « E c / K. 
As long as A > 1 the additional energy cost as com- 
pared to the state \R) is much smaller than Eq. The 
next configurations are those of a total width w c = 4 
(w c being the number of neighboring islands involved in 
the configuration), |...01 -1^010...) = |i?;l,2) and 
| ... 1 — 1r 1 . . .) = \R; 2, 1) with the charging energy 
E + E Xj2 , where E ia = U(0) - 17(1) - U{2) + C/(3) w 
2Eq I A. Thus we conclude that the regime of dominating 
charging energy Ec > Ej splits into two: 

a) Strong Coulomb blockade regime: Ec > AEj. In 
this case the charging energy difference, ~ 0(1)E C / A, 
between the charge configurations with charge-anti- 
charge pairs and the basic one | . . .001^00...) is higher 
than the tunneling energy Ej. Thus, the charge configu- 
rations of higher energy play little role. The basic charge 
configurations form a trivial tight binding band with dis- 
persion E(k) = —Ejcosk. It is this regime which was 
analyzed in 2D in Ref. fl2| . 

b) Small solitons regime: AEj > Ec > Ej. In this 
case several charge configurations hybridize with the ba- 
sic one and small solitons are formed. In what follows, we 
investigate this regime and we develop a tight binding ap- 
proach which allows us to treat this case numericall y. A 
similar approach for polarons was developed in Ref. [X 31 ] . 

To illustrate our approach we start by accounting only 
for two configurations, \R) and |i?;l,l). In Fig. [2] the 
structure of possible transitions between these states by 
tunneling of a single Cooper pair is shown. We observe 
that a tight binding situation arises again with two states 
per primitive unit cell. Instead of the cos /c-dispersion, we 
obtain the following 2x2 matrix 
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FIG. 2: Effective lattice and effective unit cell for the config- 
urations \R) and \R; 1, 1). The dashed box marks a primitive 
unit cell. The lines denote allowed tunneling between the 
configurations. 



where the second matrix accounts for the charg- 
ing energies of the states \R) and \R; 1, 1). In 
what follows we omit the common energy Eq for 

(2) 

all states. Diagonalizing Hj, yields two bands as 
shown in Fig. [3] (blue dotted curves). Next we add 
the charge states ...01 — 1^. 010...} = |i?;l,2) and 
|. . . 1 - 1 R 1 . . .) = \R; 2, 1). We find the 4x4 tight 
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In Fig.O the single particle band, the two bands of , 

and the four bands of are shown for Ec = 20i?j 
and A = 10. Here we are clearly in the strong Coulomb 
blockade regime and inclusion of the extra states only 
slightly modifies the lowest energy band. 

The idea is now to approach the intermediate regime 
AEj > Ec > Ej by extending the number of charge 
configurations. Here we went up to the total width of 
the charge configurations w c = 7 resulting in a 32 x 32 
tight binding matrix. We investigate three regimes Ec = 
10Ej, 5Ej, and 2.5Ej (A = 10). The resulting spectra 
are shown in Fig. [5] While in the strong Coulomb block- 

E(k)/Ej 




FIG. 3: (Color online) Dispersion relation for the one-state- 
approximation, i.e., a, single Cooper pair (black dashed band), 
for the two-state-approximation (blue dotted bands), and for 
the four-state-approximation (red solid bands) as described 
in the text. We chose E c = 20£j, A = 10. 



ade regime Ec > A.Ej (see Fig. [3]) the lowest band is 
very close to the — cos(fc) dispersion of a free particle, 
the shape of the lowest band in the regime of small soli- 
tons Ec < AE'j (see Fig.[4j) changes considerably. For 
Ec /Ej = 10 which is the upper boundary of the " small 
soliton" regime the lower band still has the cosine-shape 
for \k\ < it/2. For larger values of however, the band 
becomes very flat which corresponds to zero group ve- 
locity or, equivalently, to infinite mass. For smaller ra- 
tios Ec/Ej, we find that the region in the center of the 
Brioullin-zone, which is cosine-like or parabolic, becomes 
smaller (|fc| < 7r/4 for Ec/Ej = 2.5). The remaining 
flat region shows a weak oscillatory behavior. We cannot 
exclude that it is due to an insufficient number of charge 
configurations included. Indeed, while for Ec/Ej = 10 
the numerical convergence for the lowest band is good, 
it somewhat deteriorates for smaller values of Ec- For 
Ec = l.hEj the first and second bands approach each 
other at |fc| w 7r/4. This could give rise to Landau-Zencr 
transitions for an accelerated soliton. 

Soliton shape. We investigate the charge smearing in 
the regime of small solitons. For that purpose, we con- 
sider the charge-charge correlation function F(k, r—r') = 
(ipk\n r n rl \ip k ), where \ip k ) = Y,r,j a j( k ) \ R J) etkR is 




FIG. 4: Band structure for w c = 7 with parameters A = 10 
and a) E c = WEj, b) E c = 5Ej, c) E c = 2.5Ej. For clarity, 
only the 16 lower bands are shown. 
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FIG. 5: (Color online) Charge-charge correlator 
(ipk\n r n r > \ip k ) for E c = WEj, 5Ej, and 2.5Ej. In all 
plots A = 10. 



Our study was performed for infinite arrays with no 
disorder (offset charges). In the limit A » 1, both, the 
array borders and the offset charges create smooth vari- 
ations of the potential energy of a Cooper pair (wells or 
barriers). The amplitude of these variations ~ 0(1) AEc 
is, however, very large. The propagation of charge will 
thus crucially depend on the dispersion relations obtained 
in this paper as well as on the dissipation in the system. 
Further studies of these issues are necessary. 




FIG. 6: Quadrupole moment of the charge-charge correlator 
{tp k \n r n r , \tp k ) for E c = WEj, 5Ej, and 2.5Ej. For \k\ < 
7r/4 the Lorentz contraction can be observed. 

We acknowledge numerous discussions with A. Usti- 
nov, R. Schafcr, H. Rotzinger, and B. Malomed as well 
as the participants of the SCOPE 2009 meeting in Karl- 
sruhe. We thank I. Martin for pointing Ref. [13| to us. 



the Bloch wave function of the soliton (lowest band). 
Here \R,j) denotes the j-th charge configuration cen- 
tered at the island R, e.g., \R), \R;1,1), \R;1,2) etc.. 
We obtain F(k,r - r') = J2 3 l a j'(^)| 2 Cj( r ~ r ')i where 
Cj( r ) = J2 r ' n ^ n r''+r (this quantity is independent of 
R). The quantity is the number of charges on is- 
land r for the charge configuration \R,j). Note that the 
correlation function is normalized, i.e., ^2 r ,(n r n r >) = 1 
if we choose the normalization of the Bloch wave func- 
tions such that X)jl a j| 2 = 1- hi Fig-El we plot the 
charge-charge correlation function in the whole Brillouin 
zone. We observe extended structure appearing in the 
flat band regions. To characterize the width of the charge 
distribution we plot in Fig. [5] the quadrupole moment 
Q(k) = ~^2 r r 2 F(k,r). For small values of k we ob- 
serve the Lorentz contraction, as predicted by the sine- 
Gordon model. In the region of flat dispersion (infinite 
mass) the soliton becomes much wider. A question arises 
whether a model of sine-Gordon type could explain this 
phenomenon. 

Discussion. In this paper, we have identified the 
regime of small charge solitons and investigated numeri- 
cally their properties. One of the characteristic features 
is the flattening of the dispersion relation in the outer re- 
gion of the Brillouin zone and simultaneous broadening 
of the soliton. 
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